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Abstract 

We consider complex spectrum of unstable model Hamiltonian and its physically meaningful 
regions. It is shown that physically meaningful regions are obtained by taking only Hardy class 
below part of test functions. Through specific model it is shown that this complex eigenstates 
with restricted test functions show causality and represent physically meaningful regions in time 
evolution. 
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I. INTRODUCTION 



The introduction of complex eigenstate has a long history. Since Gamow proposed com- 
plex energy eigenstate to describe decay process itj_s^idely used for its simplicity and 



predicting power for decay rate. Matthews and Salam {2 4| developed the concept of an un- 
stable particle in terms of asymptotic states. Nakanishi [5] introduced complex distributions 
to define a complex eigenstate of the Hamiltonian in Lee's model joj. The contour defor- 
mation method in momenturn space is also studied and applied in nuclear physics 7H9|. 
Sudarshan, Chiu and Gorini |lO[ constructed complex eigenstates using contour deforma- 
tion in the complex plane. Bohm and Gadella ll| constructed ^complex eigenvectors using 
poles of the S matrix and Hardy class test functions (see also [12]). Prigogine and collabora- 
tors studied extensively the properties of complex spectral representations in the Friedrichs 



model 
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Moiseyev 



, and defined unstable states in Liouville space (see [l^ and references therein). 



15l | introduced complex scaling method to make diverging complex eigenfunctions 



square integrable functions. In this way the evolution of a resonance state at large distances 
is "compressed" in a small part of spaces. 

But its physical meaning has been in debate. The complex eigenstate itself is not normal- 
ized and its resemblance to the real wave function evolution ceases to exist outside certain 
spacetime region. When the decay occurs it has starting time, and the decay products 
propagate according to the causality conditions. In other words, its physically meaningful 
spacetime region is related to its initial conditions. Notable researches to this direction have 
been done by Bohm and Gadella, who proposed certain class of functions, namely Hardy 
class functions, should be applied to complex eigenstates to obtain physically meaningful re- 
sults. They also proposed rigged Hilbert spaces to distinguish functional and test functions 
[igI . This removed exponential blowup in negative time, but spatial exponential growth in 
complex eigenstate still remained. 

In this article we develop a method to find physical meanings of complex eigenstates in 
all spacetime regions. We start from a specific decay model, called Friedrichs mdoel, which 
solutions and complex eigenstates are well known. Starting from the point that complex 
eigenstates are distributions and suitable test functions are needed to give corresponding 
physical meanings, we separate test function part from complex eigenstates and time evolu- 
tion operator and look for suitable test functions which give physical meanings to complex 
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eigenstate. It is shown that in this model Hardy class functions are related to the causality 
condition and by restricting test functions to Hardy class function below we con obtain 
physically meaningful regions of complex eigenstates. 



II. MODEL 

To study the physical meaning of complex eigenstates, we consider Friedrichs model which 
complex eigenstates are well known. This model describes a discrete system interacting with 
a continuous scalar field. The Hamiltonian is given by 

/oo /'OO 
dkuk\k){k\+\ / dkvk{\l){k\ + \k){l\) (1) 
■oo J —oo 

where we put c — H—l. The state |1) represents the bare atom in its excited level with no 
field present, while the state \k) represents a bare field mode ("photon") of momentum k 
together with the atom in its ground state. In spatial representation 

m = ^e''^ (2) 



which is delta normalized. The excited state is analogous to an unstable particle state, 
while the photon is analogous to decay products. The term Xvk{\l){k\ + \k){l\) stands for 
the transition from 1 state to k state and from k state to 1 state. 

The energy of the ground state is chosen to be zero; ui is the bare energy of the excited 
level and ujk = \k\ is the photon energy. A is a small dimensionless coupling constant (A <^ 1). 
We shall consider a specific form of the interaction potential 

= i+Mur 

The constant M^^ determines the range of the interaction and gives an ultraviolet cutoff. 
Other forms of potential (form factors) may be treated in a similar way. 

Prom the dispersion relation ouk — \k\, the free- Hamiltonian eigenstates \k) and | — k) 
have the same eigenvalue ouk- We remove this degeneracy by rewriting the Hamiltonian 

/•oo /"CO 

Hf^uj,\1){1\+ dkuk{\Sk){Sk\ + \Ak){Ak\)+ dkV2Xvk{\l){Sk\ + \Sk){l\) (4) 
Jo Jo 



where 



\S,) ^ ^{\k) + I - k)), \A,) ^ ^{\k) - I - k)). (5) 



From Eq. (jl]) we see that the discrete eigenstate |1) only interacts with the symmetric field 
eigenstate l^fc). From now on, we concentrate on only the discrete atom state and the 
symmetric field states of the Hamiltonian. Rewriting 

\uj) = \Sk), v^ = V2v^^, (6) 

we get the atom-field interaction Hamiltonian 

H = Ho + XV (7) 

roQ poo 

Ho = uJi\l){l\+ dcu uj\uj){uj\, \V = \ dtu v^{\l){uj\ + \uj){l\). (8) 
Jo Jo 

It is written as 

/•oo roo 

H = ui\l){l\+ duu\uj){uj\+ dujXv^{\l){uj\ + \uj){l\) (9) 
Jo Jo 

where the first term and second term represent discrete system and continuous field each, 
and third term represents their interactions. A is coupling constants and v^j represents 
interaction for each mode u. 

This Hamiltonian has an exact diagonalized form. When ui > duX^v'^/u, We can 
write 

/»00 /"OO 

H= duu\F+){F^\= duJu\F~){F~\ (10) 
JO Jo 



where 



I \ , ^'"^ ii\ , ^'^^ Ha ' ^^^'1^') 
l^cj) = |w) + -T7-t|1) + -TT-T / 



r]+[oj) ri+[u) Jq uj — uj' + ie 



|i^j) = l-) + ^|l) + ^ frfc'^^^^ (11) 



with 



r^±(z)=z-ui- du^-^ = z - ui - P du—^±t'KX^vl (12) 
Jo -u Jo z-u 

where P J means Cauchy principal integral. We can choose + branch or — branch for the 

diagonalized solution. In Eq. f lT2|) . l/{z^ — u) means that z is analytically continued from 

above (+) or below (— ). For real z, it can be understood as = z ± ie, where e > is 

infinitesimal. These eigenstates satisfy the eigenvalue equation as well as the orthonormality 

and completeness relations 

POO 

H\F^)=u\F^), {F^\F^,) = 5iu-u'), / du^ \F^) {F^\ = 1. (13) 

Jo 



It is noted in no discrete eigenstate is present. The discrete state |1) is represented as sum 
of continuum state and it decays into continuum if there is no initial field. 

When complex eigenvalues are allowed, one can find complex eigenstates in this model. 
Explicit forms of complex eigenstates can be written as 



oo 



^.i) = <'(|l)+ / ^^^^) (14) 



'0 % 



'Z\ I 



and 



oo 



(17) 



In Eq. (HI]) - Eq. ([T7 



-oo x2 2 

iV, = ( 1 + / rfu;-p^ ) (18) 



z; — oo 



1^1 

and Zi is a complex number on the lower half plane satisfying ri'^{zi) = 0. This complex 
eigenvectors have eigenvalues 

(0, J = (0, J ;2i , I = (0,c.c 1 (20) 

If time evolution operator e~*^* is applied to complex eigenvector \zi), it gives e~*^^* factor 
which grows exponentially for negative time. Similarly time evolution operator acting on 
givesfactor which grows exponentially for positive time. 

Exponential growth also appears in spatial domain. If {x\zi) is considered, it yields 

= N^ / dco ' ^ (21) 

Jo z^ -u 

and its spatial feature is shown in figure 1. Exponential growth is shown as |x| increases. 

In this sense complex eigenstates have unphysical spacetime regions which are not shown 
in real time evolution. Next section we try to remove this unphysical features in spacetime 
by considering test function spaces. 
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III. TEST FUNCTION SPACES FOR COMPLEX EIGENSTATES 



In this section we investigate the conditions for complex eigenstate and its test functions. 
For this purpose construction process of complex eigenstate is examined first. 

In Priedrichs model, one can start complex eigenstate construction from the pole of S- 
matrix. From completeness relation is 

poo poo 

/ MF;^){F:\^ dw\F-)S{co){F:\ (22) 
Jo Jo 

where 

(23) 

This S- matrix has poles whenever r}~^{z) — 0. To extract this pole contribution, one can 
change the integration contour such that 

poo p p 

/ ch;\F-)S{u;){F^\= i d^|Fj)5H(F+| + / (24) 

^0 Jzi JC 

Here, dcu means integration over a small circle centered at Zi in clockwise direction. If 
S{uj) has a simple pole at zi, this integration gives residue value at zi. Explicitly 



/ du\Fj)S{u){F+ 

J Zl 



Zl 



r] (u) T] (u) Jq cu — u' — ie J rj'^ (u) 



V V {^) V {^) Jo UJ-U' - 16 J 

Using residue theorem and making use of the fact 

r]+{uj)^r]-{uj) + 27riX\l (26) 
Ni = hm ^ (27) 



we get 



/ duj\F-)Siuj){F:\ 

J Zl 

= /"cia;^^-27riA^;,,|^i)Vif(l|+ H dco^^^^ - 2niXvM] 

\ Jo Zi-OJ J \ Jo Zi-OJ J 

= l0.i)(0.J- (28) 
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A few points noted while deriving complex eigenstates. First complex eigensates are derived 
in dyad form, since it is derived from the completeness relation. Second the test functions 
which act on complex bra and ket should have analytically continuation from real energy 
domain. 

The physical part of complex eigenstate can be traced by considering original complete 
set. Arbitrary product {f\g) expressed by complete set 

{f\g) = r dco{f\F^){FM = r du:U\F-)S{u:){F^\g) = H du;-^ (29) 
Jo Jo Jo ^ ~ ^1 

with 

h{uj) ^{co- z^){f\F-)S{uj){FM- (30) 

Now main question is which portions of h{uj) should be contained in complex eigenstate to 
be physically meaningful. For this purpose, let us first examine the above equation closely. 
When h{u!) is square integrable and its support is in [0, oo), it can be written as 



/ duj^^ = / duj^-^ = 27rih(zi)"-. (31) 
Jo u-z J_^ u-z 

This h{zi)^~ term is called Hardy class below part of h{zi), and has useful properties. Since 
the original complete set gives Hardy class below part of h{u), it is quite tempting that we 
sue this class of functions as test functions. Before we go further it is desirable that the 
definition and properties of Hardy class functions are stated. 



A. Hardy class functions on the real line 

A complex function G(E) on the real line is a Hardy class function from above (below) 
if G(E) is the boundary value of an analytic function G{uj) in the upper (lower) half plane 
and 



J — c 



dE\G {uj)\P < oo (32) 



for all uj in the upper (lower) half plane. Here we consider p — 2 cases since in quantum 
mechanics we deal with square integrable functions. The spaces of above or below Hardy 
class functions are denoted as if? or H'^ . 



Hardy class functions has important properties that can be quite useful for our construc- 
tions. 

(1) If G{u) is in then G{uj) on the real axis is uniquely determined by its values on 
the positive real axis. 

(2) If G{uj) is in Hi then 

J. r . / «(^' ""W > ° , (33) 

2m J_^ u-z 1 for Im[z\ < 
If G(u) is in then for all Im[z] < 



(3) If is in H^, then the Fourier transform 



G±{t) = -= I dwe-"^'G{uj) (35) 



has the property 

G+(t) = for i < (36) 
G_{t) = for t > (37) 

(Paley- Wiener Theorem) 

(4) If G{(x!) is in H^, its complex conjugate is in and vice versa. 

The above properties gives quite simple way to decompose square integrable functions on 
the real line into Hardy class functions. One can write 

ffuj) r duj' r duj' ^^^'} = r dusiu - uj')fiuj') (38) 

for the positive infinitesimal e. For square integrable functions, this relation holds if f{uj) 
is continuous (even if f{oj) is not continuous discontinuous points are in measure zero set it 
might be ignorable in physical situations). When the above relation can be defined, one can 
make Hardy class decomposition on the real line as 

/(a;) = /(a;)^+ + /(c^)^- (39) 
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TABLE I: Examples of Hardy class decompostion. t > and zi is on the lower half plane. 



f{oo) 






Q—iuit 


^ — izit 


Q — iujt Q — izit 


U> — Z\ 


lO — Z\ 


UJ — Zl 









U — Zl 


UJ — Zl 



with 



2m 
1 



uj — uj' + ie' 



27rz UJ — UJ' — te 



(40) 
(41) 



Note that if e is not infinitesimal then Hardy class functions above and below can be defined 
without any problem on the upper or lower half planes. 



/(^)^^+ = 7^ duj'^^ if Im[z] > 0, 



2TTi ./_oo Z-UJ' 



f{z)^- = ^ r duj'I^ if Im[z] < 0. 
2m z-uj' 



(42) 
(43) 



The table shows examples of hardy class decomposition. One can see that Hardy class 
decomposition above (below) removes original function's divergence and singularities on the 
upper (lower) half plane. 

One can also consider Hardy class operators using Eq. ( HOj) and Eq. ( HTi) . If we define 

/(u;)^+ = [f{uj)f^ 



fiujf- = [fiuj)f- 



2Txi J_r^ UJ ~ uj' + le 



1 

2^i 



— oo 
oo 



duj 



UJ ~ UJ' — It 



(44) 
(45) 



then we have operators which take if + or H— part of test functions. They are linear and 
orthogonal if test function's if + and H— parts are defined on the real line. 



\aj{uj) + hg{uj)\^^ = aifiuj)]^^ + b[g{uj)]''^ 
[/(a;)^+]^+ = [fiuj)f^, [fiujf-f- = [fiuj)]"- 
[fiujf-]^-' = [/(u;)^+]^- = 



(46) 
(47) 
(48) 



Hardy class decompositions are frequently used in physics though its name is not often 
mentioned. When diverging denominator such as 1/ (E — H) is regularized as 1/{E — H±ie), 



it means this distribution takes Hardy class above (below) part of test functions as seen 
in Eq. f l40|) and Eq. pT]) . Calculations of retarded and advanced Green's functions and 
Lippmann-Schwinger equation use this regularization. 

The Hardy class function above (below) behaves very nicely in the upper (lower) half 
plane. It is analytic and monotonically decreases as the imaginary part of argument are 
increasing absolutely. This monotonic decreasing property can be associated with physi- 
cal properties of test functions, like considerations of boundary conditions in retarded and 
advanced Green's function derivation. Next section we apply Hardy class test function 
conditions to our model Hamiltonian and examine its physical properties. 

IV. HARDY CLASS TEST FUNCTIONS AND COMPLEX EIGENSTATES 

In this section we apply Hardy class functions as complex eigenstate dyad and observe 
its physical properties. In Friedrichs model the complex eigenstate dyad is given by 

\K){<P.. I = iVi (|i) + [ rf-l^^) + [ ^-|^) • (49) 

Applying (/|, e~*^* and \g) on complex eigenstate dyad, one gets 

= A^if(/|l)+ r du^^^^^^-27^^\v,M\zA 

V Zi-UJ J 

xe--*((l|^7) + /"^^^^^ - 27rat;,,(zi|^)). (50) 

In Eq. fl50|) (/l-Zi), (-^ilfi') and e~*^^* give exponential divergences in spatial and time domain. 
To obtain physically meaningful regions these terms should be modified. Also it is desir- 
able that complex eigenstate is not modified while obtaining physically meaningful results, 
otherwise we mix states with initial conditions. 

Mathematically choosing Hardy class part is quite promising. Exponential divergence 
happens when field part of test functions are analytically continued from real u to the 
complex Zi in lower half plane. If Hardy class below parts are taken, they no longer diverges 
since H— parts monotonically decreases as its imaginary part decreases and analytic in the 
lower half plane. 

By expanding Eq. ( l50l) and taking H— parts only on test function part (thus keeping the 
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form of complex eigenstates), the modified {f\e~ 



iHtl 



\g) can be written as 



- modified 



+N^{-27r\v,,)[{f\Lo)e 



Zi — uo 

'iu)t~\ H— 
Izi 



ir (ii9> + 







'Mg) 



,Xv^:{u'\g) 



) Zi — w 

+Nr{-27:\v,,f[{f\uj)e''^\uj 



duo 

2:1-0;' 



(51) 



where means Hardy class below part of F{(jj) at w = ^i. 

In Eq. (!5T]) H— parts at Zi are taken for four different terms; e"""^*, (/|u;)e^*'^*, e^*'^*((X'|(?) 
and (/|(X')e^*'^*((X'|(?). Close examination reveals quite interesting physical meaning to Hardy 
class below parts. 

[e"*"^*]^" is related to |1)(1| which is projection operator to the discrete state. It yields 



-iu}tT^H- 



0(t) e 



-iz\t 



(52) 



where 



m = { 



1 for t > 
i fort = 
for t < 0. 



(53) 



This shows exponential decay only for positive time (for t = test function is constant and 
not normalizable) . Exponential growth in negative time is removed and only decay to the 
future remains. 

[(/|a;)e~*'^*]^~ is related to which is discrete to decaying field transition. Its 

physical meaning becomes clear when (/| is chosen as (x| with x 7^ 0. For a; = case only 



-iujt 



term remains and the result is the same as previous one. It yields 



[{x\u)e' 



-iLjf\H- 
izi 



1 



2v^ 



(0(t - |x|) e-'"^(*-'^') + Q{t + \x\) e-^^i(*+l^l)) 



(54) 



The result shows causality of decaying field. Only parts of incoming field e-*^i(*+l^l) which 
can affect |1) for positive t are selected, which are in the region t < —\x\. For the outgoing 
field e~*^i(*~l^l) only the parts which are produced by discrete state for positive t are selected. 
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[e~'^'^^ {uj\g)]^~ is related to which is decaying field to discrete state transition. 

When l^f) is chosen as |x) it yields the same result as [(x|a;)e~*'^*]^~ since {x\u) is real. 
[{f\uj)e~'^'^'^{uj\g)]^~ is related to field to field transition. When (/| and \g) are chosen as 
and 1 0:2), it yields 

^{Q{t+\xi\ + |x2|)e-'^i(*+'"^'+l"2l) + e(t- |a;i| + |x2|) e-'^^(*-l"il+l"^l) 

+e(t + |xi| - |X2|) e-^^i(*+l"il-l"2|) + -\xi\- \X2\) e-^^i(*-l-il-l^'^l)^ . (55) 

Causality relations are also shown in Eq. fl55|) . Without loss of generality we can assume 
|xi| > \x2\- e"*^^*^*^'^^'^'^^'-* shows decaying incoming wave for Its support is t > 

— |xi| — |x2|, and only parts of incoming field which can reach from —\x2\ to (or from 
\x2\ to — |a;i|) for t > are selected. e-*^i(*-l^il+l^2|) g^ows decaying outgoing field for 
Its support is t > |xi| — IX2I, which is the time needed for a field at \x2\ and t > travels to 
On other hand, e"*^^^*^!^'^!"!^^!) is incoming wave for Its support is t > — |xi| + |x2|, 
which is the time that incoming wave at Xi travels to \x2\ for t > so that transition from 
|xi| to IX2I is made, e"*^!*^*"!^'^!"!^'^!) is outgoing wave for |a;i|. Its support is t > + \x2\^ 
which is the time needed for a field at \x2\ and t > travels to xi. 

All above relations show causality for both discrete state and field in regard to t = 0. 
The results show that choice of Hardy class functions as test functions of complex eigenstate 
selects only causal part of complex eigenstate and removes unphysical divergences. 

V. CONCLUSION 

We studied complex spectral representation and it test functions in Friedrichs model. 
The choice of Hardy class functions as test functions of complex eigenstate restores causality 
relations and removes divergences of complex eigenstate in spacetime domain. 
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